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ABSTRACT 

Several neutral hydrogen (HI) cavities have been detected in the Milky Way and other nearby star 
forming galaxies. It has been suggested that at least a fraction of them may be expanding supcrshells 
driven by the combined mechanical feedback from multiple supcrnovae occurring in an OB association. 
Yet most extragalactic HI holes have neither a demonstrated expansion velocity, nor an identified OB 
association inside them. In this work, we report on the discovery of an unbroken expanding HI 
supershell in the nearby spiral galaxy M101, with an UV emitting OB association inside it. We 
measure its size (500 pc) and expansion velocity (20 km/s) by identifying both its approaching and 
receding components in the position- velocity space, using 21 cm emission spectroscopy. This provides 
us with an ideal system to test the theory of supershells driven by the mechanical feedback from 
multiple supernovae. The UV emission of the cluster inside the supershell is compared with simulated 
spectral energy distribution of synthetic clusters of the appropriate age (~ 15 Myr). The observed 
UV flux is found to be consistent with an association of the appropriate mass (~ 10 5 M Q ) and age 
required by the energy budget of the supershell. Properties of this supershell and another previously 
reported in the same galaxy are used to infer its neutral hydrogen scale height and mean neutral 
hydrogen density in the disk. The presence of another UV emitting stellar association in over-dense 
swept up gas is discussed in the context of propagating star formation. 

Subject headings: supernovae: general — Galaxies: individual (M101) — ISM: Bubbles — Radio lines: 
ISM 



1. INTRODUCTION 

OB associations and other young stellar clusters con- 
taining a few to thousands of massive stars of spectral 
class O and B; harbor the majority of massive stars in 
spiral galaxies. Stars of spectral class O are short-lived, 
and end their lives as supernovae within a few millions of 
years. These associations last for few tens of millions of 
years because all their massive stars exhaust their nuclear 
fuel within this time. Energetic radiation from OB asso- 
ciations ionize the interstellar medium (ISM) in galaxies 
and produce localized HII regions as well as diffuse ion- 
ized gas. Supcrnovae of various types explode from the 
massive stars that form in these associations wh ich in 
turn pressurize and churn the ISM (|Oev et al.ll2002f ) . Ex- 
panding neutral hydrogen structures , called super shells 
were identified in the Milky Way by iHeilesI (| 19790 . Su- 
pershells form and evolve due to the continuous mechan- 
ical energy injection by s tellar winds and SNe in m assive 
parent OB associations (jMcCrav fc Kafatoslll987t ). The 
outer shocks of the superbubbles sweep up the ambient 
ISM into a thin cool shell. 

In this work we report yet another striking case of 
an expanding HI supershell in M101 which we argue 
is driven by multiple supernovae from a young stellar 
association. We use the observed properties of the HI 
supershell to determine its dynamical age and energy 
requirement. Here we demonstrate that this HI cavity 
has measured expansion velocity and at the same time 
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is fully contained, i.e. it has not yet broken out of the 
disk of its host galaxy. We also detect a UV emitting 
stellar association inside the HI hole. We model its spec- 
tral ^ner^v_dis_tribjrtfon ; _te^ and confirm the hypothesis 
of (jMcCrav fc Kafatosl fl987h that some supershells are 
driven by the combined action of multiple supernovae 
occurring in a stellar association inside it. 

These supershells may also play a role in the process of 
star formation and favor propagati n g star formation, as 
pointe d out by IMcCray fc Kafatosl (|1987D : IPalous et all 
(|1994f K iHopkins et all (|2008| ). point out that the neutral 
and molecular gas replenishment in the walls of the su- 
pershells may provide the trigg er for collapse and fur ther 
star formation. In NGC 300, iBlair fc LonH (|1997t) re- 
ported several cavities, much larger than g alactic super- 
nova remnants (SNRs). iPayne et al.1 ()2004l ) suggest that 
these may represent the above mentioned supershells. 

Cavities have been found in neutral hydrogen (HI) 
surface density maps of several nearby galaxies (for 
a classical review and inventory of HI holes, sec 
iTenorio-Tagle fc Bodenheimen fl988l for a recent review 
and survey of HI holes, see IBagetakqs et aT] 120101) : in- 
cluding M101 ([Kamphuis et all ll99lf T However, the 
presence of a cavity alone is insufficient evidence for 
the existence of a supershell. It may simply be a low 
density region between higher density ones, created by 
turbulent distribution of gas. This pitfall is illustrated 
by the case of Holmberg II. In that galaxy IPuche et al.l 
(|1992h used underdensities in HI column density map 
to create a catalog of HI cavities and postulated that 
they have b e en ev acuated by multiple SNe. However 
iRhode et al.l (|1999[ ) rule out a multiple supernova origin 
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for most of these holes since they were unable to find 
any young cluster or even trace of recent star formation 
activity which would have led to the SNe responsible for 
these Hl-hole s. iStewart et al.l (120001) have since then re- 
examined the lRhode et al.1 (|1999h results and concluded 
that the energy available from massive stars inside HI 
shells indicates that energy deposited into the ISM from 
supernovae and stell ar winds is s ufficie nt to account for 
the HI morphology. iSilich et al.l (|2006() have found that 
the observed properties of the most prominent kpc-scale 
structure in IC 1613 and the level of the detected star 
formation activity are inconsistent with the hypothesis 
that they are formed by energy injection from multiple 
supernovae. 

The lack of direct determination of the source of energy 
driving these bubbles in many case s have led to numerous 
suggestions regarding their origin. lLoeb fc Pernal ffl998h 
have suggested that a fraction of HI super shells might 
be dri ven by Gamma Ray Bursts (GRBs). IWada et al.l 
(|2000() have discussed the possibility of forming cavi- 
ties from the nonlinear growth o f thermal and gravi - 
tatio nal ins t abilit ies in the ISM. iTenorio-Taglel (| 19811 ) 
and iHeilesi (|1984l ) have suggested that the impact of 
small companion galaxies or high velocity clouds with 
the disk of a larger galax y may form these cavities. 
lElmegreen fc Chiang] ( 1982f ) have proposed that radia- 
tion pressure from field stars may p rovide a driving mech- 
anism for supershells. lOevI ()1996[ ) has provided a com- 
prehensive comparison of the supernova driven supershell 
model and the observed properties of the LMC shells . 
Outstanding issues in the field are listed byEiyl (l200l . 

The first unambiguous case of an extragalactic expand- 
ing HI shell associated with a hole in the HI su rface den- 
sity was reported by iKamphuis et al.1 ([l991) in obser- 
vations of M101 using the W esterbork Synthesis Radio 
telescope. iWeisz et al.l (|2009| ) have investigated the cre- 
ation of a supergiant HI shell in a dwarf irregular galaxy 
IC 2574 and its role in triggering star formation around 
it. Star formation studies in most nearby galaxies in- 
dicate that sites of recent star formati on are correlate d 
with sites of higher HI surface densities (jKennicu tt 1998), 
leading to several recipes for star formation. The feed- 
back processes described above should however clear out 
the gas and lead to young stellar associations in regions 
of low gas density surrounded by higher density swept 
up shells. Therefore we began searching for cavities in 
the HI surface density, occupied by young clusters. 

2. DATA DESCRIPTION 

In this work we have used 21 cm line data to identify 
and measure kinematic properties of our putative super- 
shell. We have further used UV data to identify the 
young stellar association responsible for driving it. We 
discuss below, the detailed nature of the data; especially 
the sensitivities, spatial and spectral resolutions of the 
two relevant surveys. 

2.1. The HI Nearby Galaxy Survey 

The HI Nearby Galaxy Survey (THINGS) undertaken 
at the NRAO Very Large Arra y (VLA) to study t he HI 
emission from nearby galaxies dWalter et al.l [2008h pro- 
vides high angular resolution (~ 7") maps of HI emission 
and high spectral resolution (~ 5km/s) data. The sur- 
vey has acquired radio interfcromctric observations for 



the target galaxies and presents the data in the form of 
image cubes generated from the interfcromctric data af- 
ter background continuum subtraction. The three axes 
in these cubes represent the Right Ascension (RA) , Dec- 
lination (Dec) and recession velocity. For each triad of 
these coordinates, the cube gives the HI line flux den- 
sity, in excess of the best estimate for the continuum. 
Information about the position and velocity of neutral 
hydrogen gas can be extracted from this data cube. 

Demanding that the diameter of the supershell be more 
than the projected beam size, we can only detect super- 
shells with radii more than 24 pc (#io-Di), where #io is 
the beam size in units of 10" and D\ is the distance to the 
host galaxy in Mpc. The two velocity components can be 
resolved only if the velocity is greater than the spectral 
resolution. This is of course limited by the sensitivity to 
HI in the beam, which is about ~ 10 20 atoms/cm 2 for 
the data being used. 

2.2. GALEX Nearby Galaxies Survey 

Data on M101 from Far and Near Ultraviolet ob- 
servations are available from GALEX Nearby Galaxies 
Survey, co nducted with the G alaxy Evolution Explorer 
(GALEX) (jMartin et al.ll2005[ ). The ultraviolet emission 
provides us with a simple tool for identifying young stel- 
lar associations, which are expected to host a sizable frac- 
tion of the core collapse SNe, which supply the energy 
budget for the expanding supershells. We use informa- 
tion about both the HI distribution and the young stellar 
populations, to search for neutral hydrogen supershells 
powered by multiple supernovae. 

The GALEX data is provided as FITS files where each 
pixel gives the counts per second (cps), in the relevant 
band, for each RA-Dcc pair. Comparing the HI sur- 
face density maps with those for near ultraviolet (NUV) 
emission for M101, reveals multiple regions with high 
UV emission but low HI surface density, surrounded by 
regions of higher HI density. These include the HI su- 
perbubble reported by Kamphuis et al.(1991). Another 
association hosting SN 1970G, is detected in the ul- 
traviolet and found to be located on the high column 
density HI ring surrounding an undcrdensity. We have 
also found a separate cavity in the HI surface density, 
hosting a young stellar association of NUV magnitude 
rriAB = 17.37 ± 0.08, surrounded almost completely by 
higher HI surface density regions (Fig. 1). This region, 
which shall be discussed further in this work, is located 
approximately at 14h 02m 25s +54d 27' 45" . 

3. EVIDENCE FOR EXPANSION 

Most extragalactic HI cavities have no demonstrated 
expansion velocities (see exceptions in Introduction) 
Demonstration of expansion velocities is essential for the 
identification of HI cavities associated with supershells. 
Even where expansion is claimed, the measured velocity 
is sometimes close to the typical rms velocity. In absence 
of clear detections for both forward and backward veloc- 
ities, claims for ex pansion are ambiguous. For example 
IWeisz et al.l <|2009f ) infer a velocity from a break in the 
P-V diagram, which they claim to be an indirect mea- 
sure of expansion, yet consistent with presently stalled 
expansion. Because of sensitivity limitations, the veloc- 
ity profile of the HI g as in an unde r dense region becomes 
difficult to measure. iPuche et al.l (|1992f) point out that 
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because of the low brightness contrast between the shells 
and the face of their Hl-holes, they find it very difficult 
to look for expansion in position velocity (P-V) diagrams 
and to demonstrate both incoming and receding compo- 
nents. The high sensitivity of the THINGS survey, used 
in our study, overcomes this limitation. 

To study the dynamics of the HI gas, we extract 
the relevant region (entire area covered in Fig. 1) 
out of the THINGS data cube and transform it, us- 
ing the Astronomical Imag e Processing System (AIPS 
|http : / / www. aips .nrao.e du/), to obtain velocity informa- 
tion on the sky plane. We slice the resulting data sub- 
cube along nine equally spaced Dec values over the entire 
region of interest (whole of Fig. 1). The HI line intensity 
is then plotted in the P-V diagrams for these slices (Fig. 
2) The P-V diagrams have a spatial resolution of 10" 
and a velocity resolution of 5 km/s. The HI distribution 
in the starting and ending slices (Fig. 2) are continu- 
ous over the RA and velocity, indicating that most of 
the gas is at the local rest frame velocity dictated by the 
galaxy rotation curve. The Dec slices cutting across the 
HI underdensity however indicate a clear transition to 
a two-component HI emission. The true nature of the 
Hi-hole can however only be studied from its detailed 
velocity map. The data slice along Dec slice D clearly 
shows that the HI under-dense region is characterized by 
reduced emission from the rest frame velocity of the HI 
disk. This provides an unambiguous detection of an ex- 
panding HI supcrshcll surrounding a UV emitting stellar 
association. The difference in velocity of the two com- 
ponents in the position velocity diagram (Fig. 2 panel 
D) is around 40 km/s. Interpreting this as twice the ex- 
pansion velocity of the supershell, we get Vs ~ 20 km/s. 
The diameter of the cavity, which lacks HI emission at 
the rest frame velocity of the disk, is around ~ 30", cor- 
responding to Rs ~ 5 00 pc at th e dista nce of 6.85 Mpc 
to M101 reported by iSaha et ail (|2006[ ). We also note 
the bilateral symmetry in the emission from the two de- 
tected components and suggest that the supershell must 
have started out at or near the mid-plane of the HI disk 
of the host galaxy. 

We also overlay the HI column density contours on a 
map of rms deviations of velocities from the mean, ob- 
tained by taking the second moment of the velocity distri- 
bution of the emission at each pixel (Fig. 3). The square 
of the rms velocities trace the line of sight kinetic energy 
per unit mass of the HI gas. Hence the high rms veloc- 
ity of up to 18.1 km/s in the Hi-hole region is evidence 
for recent energy injection into the ISM. For comparison, 
the neighborhood of the Hi-hole, reveals an average rms 
velocity of 10.5 km/s with a standard deviation of 1.4 
km/s. We construct a continuum subtracted HI spec- 
tra (Fig. 4), integrated over a box (2.6xbeamsize on 
each side) inside the Hl-holc. It clearly shows that there 
are two components of HI emission with easily differen- 
tiated velocities. We interpret the two components, as 
demonstrated in Fig. 2 and Fig. 4, as the approaching 
and receding hemispheres of an expanding HI supershell 
around an evacuated bubble. Note that the sensitivity 
of the THINGS survey allows us to detect the gas even 
in the lowest density region with high significance. 

4. IMPLICATIONS FOR HOST GALAXY 



HI observations of facc-on galaxies map their neutral 
hydrogen surface densities. However HI scale heights 
have been directly measured mostly for edge on systems. 
Scale height measurements for face-on galaxies are rare. 
Here we propose to put constrains on the HI scale height 
of M101 based on the propertie s of the described super - 
shell and the one repo rted by iKamphuis et al.l (|1991|) . 
iBruhweiler et al.l (|1980l ) have discussed the sizes of galac- 
tic HI shells as a function of the neutral hydrogen scale 
height. 

4.1. Supershells in Ml 01 

IKamphuis et al.l (|1991[ ) have argued on the basis of the 
size of the hole found by them and the velocity width of 
its HI profile that it has already reached a state where 
parts of it are accelerating away from the disk and into 
the halo. At the quoted diameter of 1.5 kpc, if the super- 
shell has burst out of the galactic HI disk, the HI scale 
height is likely to be less than Rs ~ 750 pc. On the other 
hand both the approaching and receding components of 
the HI profile for the hole reported in our work (Fig. 2 
panel D) are continuous in positio n and velocity (unlike 
Fig. 3 in IKamphuis et all (|199lD ) and appears not to 
have burst out of the HI disk. Therefore the scale height 
is likely to be more than Rs ~ 500 pc. 

4.2. Scale Height of M 101 

iMcCrav k, Kafatoi (|1987|) explain that the pressure 
driven phase of a supershell terminates when Rs be- 
comes comparable to the HI scale height. Thereafter, 
if the velocity of the shell at one scale height is greater 
than the typical RMS velocity of the HI, portions of 
the shell traveling perpendicular to the disc accelerate 
into the halo, fragment due to Raylcigh- Taylor instabil- 
ities and discharge the internal pressure into the galac- 
tic corona. The portion moving along the disc, enters 
the pressure-less snowplow phase and continues to evolve 
slowly (Rs oc i 1 / 4 ). Hence, the radii of the largest ob- 
served Hl-holes in face-on disk galaxies are expected to 
be slightly larger than the H I scale height. 

Observations of these two (iKamph uis et al.l (|199lD and 
this work) supershells constrain the HI scale height of 
M101 to be between 500 pc and 750 pc. The average 
line integrated flux, from the region covered in Fig. 1, 
is 0.18 Jy/beam x km /s. We use the conversion from 
IRohlfs fc Wilson! (1200(1 to get an average HI column den- 
sity of Nhi ~ 1.8 xlO 21 atoms/cm 2 . Even the least dense 
area in Fig. 1 has Nhi ~ 0.35 x 10 21 which is well above 
the detection threshold. Using a scale height of 600 pc 
consistent with the above constraints we infer the initial 
gas density by dividing the column density by twice the 
scale height to get no ~ 0.5 in units of atoms cm~ 3 . We 
have used this value while estimating the energy budget 
of the supershell later in this work. 

5. INVERSION OF MCCRAY & KAFATOS MODEL 

Massive stars in a young stellar association end their 
lives as core collapse supernovae. In an OB association 
this may result in the formation of an expanding neu- 
tral hydrogen supcrshcll driven by the mechanical feed- 
back from the com b ined action of multiple supernovae. 
IMcCrav fc Kafatosl (|1987f ) demonstrate that the rate of 
SNe from a cluster of typical IMF, varies only slightly 
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over the first 50 Myr and can be thought of as a continu- 
ous energy injection over such a timescale, powering the 
expansion of the supcrshcll. This provides a particularly 
simple framework, which we adopt, for understanding 
the evolution of a supcrshcll, pressure-driven by feedback 
from multiple supernova. The radius (i?s) and velocity 
(Vs) of an expanding supershell, given by them, are 

R s = 97 pc (iV*£;. 5 i/™o) 1/5 i 3 / 5 



and 



(N*E 51 /n ) ~ 1.1 x 10 3 . 



(6) 



and 



V s = 5.7 km s- 1 (N*E sl /n ) 1/5 t 7 2/5 . 



(1) 
(2) 

Where, 2V* is the number of massive (> 7M ) stars in the 
association responsible for creating the supershell. E 5 i is 
the energy produced by each supernova explosion, tiq is 
the assumed uniform atomic density prior to the forma- 
tion of the supershell and t? is the age of the supershell 
in units of 10 7 years. For an alternate formulation of 
the problem, which directly relates the evolution to the 
supershell to the mechanical luminosity of the driving 

cluster, see Appendix. 

The (jMcCrav fc Kafatosl fl987l ) model is invertible in 
principle, because the Rs and the Vs are directly ob- 
servable from the HI data. In this work, we therefore 
re-frame the equations as 



and 



*7 = (i?s/97pc)(y s /5.7kms" 



(N*E 51 /n ) = (Rs/97 pc) 2 (Ks/5.7 km s 



-1\3 



(3) 



(4) 



to express the variables, which characterize the super- 
shell, purely in terms of the observable quantities. These 
equations will allow us to estimate the age and mass of an 
association which will suffice to reproduce the observed 
size and expansion velocity of any given supershell. If 
there is no young stellar association interior to an ob- 
served expanding neutral hydrogen structure or an as- 
sociation smaller th an what is demanded by these set of 
equations, then the ([McCrav fc Kafatoill987t ) model can 
be ruled out. If on the other hand, one finds such a suffi- 
ciently luminous association, the model can then explain 
the origin of such a supershell. 

6. AGE AND ENERGY BUDGET OF THE M101 
SUPERSHELL 



iWeisz et al.l (|2009() have used the model of iChevalierl 
(|1974[) . to compute the energy required to evacuate the 
HI mass inside the supershell that they study in IC 2574. 
Chevalier's model deposits the total energy at the start 
of the simulation and he notes that this soon relaxes to 
an approximation of the Sedov solution. This was in- 
tended for simulating a single SNR. However this is not 
a suitable description of the energy input from multiple 
SNe as stars of different initial masses, explode at differ- 
ent points in time. In addition note that the power law 
index of time evolution (3/5 in Eq.l), due to continuous 
energy i njection, is di f ferent from the Sedov-like solution 
(used bv lWeisz et al.l (120091)') a nd is similar to the wind 
solution of lWeaver et al.l (|1977l) . 

The appropriate evolution of the supershell can be ob- 
tained by the self consistent solution of Eqs (1) and (2). 
The observed size and expansion velocity of the super- 
shell will be satisfied in this manner for 



t 7 ~ 1.5 



(5) 



Assuming E$i ~ 1 and uq ~ 0.5 (as already demon- 
strated), we need the driving cluster to be massive 
enough to harbor ~ 550 supernova yielding massive 
stars. Assuming the total mass of stars required to be 
formed to pr oduce each core collapse s upernova as 196 
M Q following iMcKee fc Williami (fl997h . we need a clus- 
ter of initial mass ~ 10 5 M Q which has been driving this 
supershell into the ISM for the last ~ 15 Myr. Only 
~ 15% of the energy from the already exploded super- 
novae (at tf = 1.5), needs to be in the kinetic energy 
of the supershell, assuming np = 0.5. This is close to 
the estimate from iMcCrav fc Kafatosl (119871 ). Accord- 
ing to Starburst99 (jLeitherer et all[l999l ) the number of 
ionizing photons with wavelength below 912A, drops to 
4 x 10 -4 times the initial rate, by the time a typical asso- 
ciation reaches 15 Myr. Hence, given the estimated age 
of the supershell it is no surprise that it is made up of 
neutral hydrogen. 

7. UV DIAGNOSTIC FOR MASS OF THE CLUSTER 

Most sup ershells in the Milky Way suffer fro m large 
extinctions (|Tenorio-Tagle fc Bodenheimerlll988l ). Even 
for many extragalactic holes, yo ung stellar popula - 
tions have not been detected inside (|Rhode et al.lll999l) . 
Us ing Starburst99 with the parameters as described 
in iLeitherer et al.l (|1999| ) we simulate the spectral en- 
ergy distribution (SED) of an instantaneous starburst 
(Salpeter IMF with a = 2.35 between 1 - 1OOM ) of so- 
lar metallicity [Z = 0.02) aged 15, 50 and 200 Myr, over 
a range of wavelengths from 500A to 3500A (see Fig. 5). 
The effective collecting area of the GALEX NUV and 
FUV arc plotted on the same graph for comparison. At 
the effective wavelength of the GALEX NUV band the 
luminosity of the 15 Myr (chosen to match the dynamical 
age of the supershell) starburst is given as 1.17 x 10 32 ergs 
sec -1 A -1 Mg 1 . It is to be noted, that as an association 
ages, over the time range represented in the figure, the 
flux in both GALEX bands falls by more than one order 
of magnitude. This motivates the use of GALEX images 
in locating young stellar association in nearby galaxies, 
which may be the source o f multiple supernovae drivi ng 
supershells as suggested bv lMcCrav fc Kafatosl (|1987l) . 

In this work, we measure the NUV emission of the clus- 
ter interior to the HI cavity using archival GALEX obser- 
vation. We find the NUV magnitude, uncorrected for ex- 
tinction, to be niAB = 17.37±0.08 . Using only th e galac - 
tic extinction of A v = 0.028 from iSchlegel et all (| 19981) . 
and scaling it to th e NUV band using the relationship 
from iCardelli et all (|1989f ). we get A NUV = 0.083, for 
a typical Milky Way ratio of total to selective absorp- 
tion of Ry = 3.1. At the adopted distance to M101, 
this gives an NUV luminosity of 1.51 x 10 37 ergs sec -1 
A -1 . Comparing the simulated SED from Starburst99 
(as above) and the observed GALEX flux, we infer the 
mass of the putative young stellar association respon- 
sible for evacuating the observed supershell. Assuming 
that the age of association is same as the dynamic age of 
the supershell as determined in this work, the association 
had an original mass of ~ 1.3 x 10 5 M Q . This is consis- 
tent with the minimum mass requirement of ~ 10 5 A/q 
that was derived from the energy requirement in the 
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iMcCrav fe Kafatosl (| 198 71) model. The inferred mass of 
this cluster is consistent with masses of the most massive 
of young stellar cl usters (or super s tar clusters) seen in 
the local universe flWhitmorel 12003ft. Thus we find that 
the IMcCrav fe Kafatosl ()1987| ) model provides a simple 
and sufficient explanation of the observed supershell in 
M101. 

8. INSTABILITIES OF THE EXPANDING SUPERSHELL 

Apart from the NUV emission from stars inside the 
bubble, we also note the even stronger NUV emission 
from the high density swept up ring to the south of the 
cavity (Fig. 6) at around 14h 02m 25s +54d 27' 15". 
We suggest that this is evidence for induced star for- 
mation at the edges of expanding supcrshell s driven by 
multiple sup ernovae from an OB association. lElmegreenl 
(|1989l[l994l ) has investigated the gravitational collapse of 
decelerating shocked layers in OB associations and pro- 
posed a condition for pr opagating star fo rmation due to 
large expanding shells. iVishniaa (|1994l ) has presented 
hydrodynamical stability analysis for cold gas bounded 
b y shocks. 

lEhlerova et al.l (|1997l) performed a linear analysis of 
the Elmegreen- Vishniac instability and derived an in- 
stability time (tb) which marks the onset of gravita- 
tional instability of a self gravitating expanding super- 
shell. Re-casting t heir expression in the notation of 
IMcCrav fc Kafatosl (fl987h . 



t b = 47 Myr (N t E 51 } 



-1/8 5/8 -1/2 -1/2 
c l n A*1.3 



(7) 



where c\ is the sound speed in the shell in units of 1 km/s 
and /xi.3 is the mean molecular weight in units of 1.3 amu. 
Using TV* ~ 5.5 x 10 2 , no ~ 0.5 and unity for the rest of 
the dimensionless quantities, we get tb ~ 30 Myr. Com- 
paring this with the estimated age of the supershell (Eq 
5), we find that the supershell is not yet prone to gravita- 
tional instabilities leading to star formation. Since there 
is star formation at the boundary of the region discussed 
before, this then may be due to the expanding supershell 
encount ering an existing d ense cloud triggering star for- 
mation (jElmegreenl (|1998| ) and references therein). Such 
a scenario can also provide a natural explanation for the 
deviation from spherical symmetry in the Hl-holc. 

9. DISCUSSIONS 

We have reported the discovery of an expanding HI su- 
pershell in M101, which hosts a young stellar association 
in its cavity. The supershell may have been created by 
the combined effect of multiple supernovae. This allowed 
us to use the McCray & Kafatos model to estimate the 
age and mechanical luminosity required from the clus- 
ter that harbored the massive stars responsible for these 
explosions. By comparing the observed GALEX NUV 
emission of the cluster found inside the Hi-hole, with the 
synthetic SED of the required age obtained using Star- 
burst99, we demonstrate that the cluster is heavy enough 
to host enough massive stars which explode as super- 
novae and provide energy to drive the supershell. We 
have discussed the implications of the measured proper- 
ties of the supershell in the host galaxy to estimate the 
galaxy's HI scale height and the mean neutral hydrogen 
density in the disk. We also note the presence of young 
stellar associations near the edge of the observed super- 



shell and discussed the role of instabilities in propagating 
star formation. 

Further studies, e.g. high resolution imaging or inte- 
gral field optical spectroscopy of the cluster in the HI- 
hole, may determine the metallicity enrichment of the 
gas and through the use of spectral fitting of stellar pop- 
ulations, allow the determination of true age and initial 
mass of the cluster. Spatially resolved integral field spec- 
troscopy also has the potential to reveal the role of multi- 
ple stellar populations, if any, within the cavity. We have 
shown that the observed UV fluxes of the cluster are con- 
sistent with the standard model of supcrshcll formation. 
An independent determination of the age of the cluster 
will prove the applicability of the model in this case. Sen- 
sitive X-ray observations may reveal emissio n from pos- 
sible hot c oronal gas inside the s upers hell ([Silich et al.l 
[2WII2005T ). However. IChuet al.l (fl995T) suggest that an 
intermittent process such as hidden supernova remnants 
may be responsible for the X-ray-bright super bubbles 
and hence, not all superbubbles are X-ray bright. 

In case of triggered star formation, iStewart et al.l 
(|2000f ) have emphasized the importance of establish- 
ing that stars are formed in chains of progressing age. 
We predict that spectroscopy will reveal the populations 
noted near the edge of the supershell, to be younger than 
the cluster seen inside the supershell. A detailed study of 
this supershell with the upcoming ALMA may render CO 
line maps which will trace the c old molecular gas at the 
sites of triggered star formation. Bagcta kos et al.l (|2010f ) 
have now produced a catalog of more than 1000 HI holes 
in 20 nearby galaxies from the THINGS sample. A lot 
of these galaxies have deep GALEX observations, which 
we have demonstrated as a simple tool for identifying 
possible young stellar associations supplying the energy 
requirement. This could form the basis of a much larger 
study to determi ne the fraction of HI hole s that can be 
explained by the IMcCrav fe Kafatosl (|1987t ) model. 
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Fig. 1. — HI surface density map of M101 at the location of 
the cavity surrounded by higher surface density of gas. Note the 
beam size is 10", approximately the same as the central under- 
dense contour. Alphabets beside the declination markings, denote 
the positions of the declination slices in Fig. 2. 



Fig. 3. — HI surface density contours, plotted over the 2nd mo- 
ment of the HI velocity distribution map (RMS deviations in veloc- 
ity from the mean) in color. Note that the highest kinetic energy 
per unit mass are observed within the Hl-hole. 



P-V diagrams along 9 DEC slices 
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Fig. 2. — HI gas velocity distribution vs RA in nine equally spaced 
Dec slices across the region shown in Fig. 1. Note the two resolved 
velocity components surrounding the Hl-holc, e.g. in panels C and 
D. 



Continuum subtracted HI spectrum 
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Fig. 4. — Continuum subtracted HI spectrum integrated over 
a box, 2.6xbeamsize on each side, inside the Hl-hole. Note the 
double peaked emission profile. 
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Fig. 5. — Spectral Energy Distribution (SED) of a 10 6 Af Q as- 
sociation (with Salpctcr IMF and solar metallicit y) at the age of 
15, 5 and 200 Myr computed using Starburst99 l ILeitherer et al,l 
[1999T) . plotted for comparison with the effective collecting area of 
GALEX satellite telescope in the NUV and FUV bands. Note that 
the flux in both bands falls by more than an order of magnitude in 
the depicted range of ages. 



COLORS: GALEX NUV Count Rate 

CONTOURS: HI Column Density 

100 200 




14 02 29 28 27 26 25 24 23 22 21 20 
RIGHT ASCENSION (J2000) 
Color scale flux range= 0.0 280.0 Milli CPS 
Contour peak flux = 5.1 033E+02 JY/B*M/S 
Contours = 5.000E+01 * (1, 1.414, 2, 2.828, 4, 5.657, 8) 



Fig. 6. — HI surface density contours, plotted over the GALEX 
NUV photon count rate (in color). GALEX typically produces 
images with 5" FWHM. Note the cluster inside the HI- hole and 
also the brighter one on the swept up shell in the south. 
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APPENDIX 

A. SUPERSHELL DYNAMICS: ALTERNATE FORMULATION 

In this work we have used the early formulation of the supershell dynamics as given bv lMcCrav fc Kafatosl (fl~98l . 
This formulation depends directly on terms like N^E^i whic h relate the dynamics to th e total number of supernova 
yieldin g stars in the associations. However in a later work, iMac Low k, McCravl (|T988) ; Bi snovatvi-Kogan fc Silichl 
(|1995D have directly linked the evolution o f the supershell with the co mbined mechanical luminosity of the stellar 
winds and supernovae from the association. IMac Low fc McCravl (|1988l ) give the radius and velocity evolution of the 
supershell as 

i? s - 267 P c(i 3 8/«o) 1/5 i? (Al) 

and 

V s = 15.7 km s- 1 (L 38 /« ) 1/5 if V \ (A2) 

where the new term is the m echanical luminosity of the cluster in units of 10 38 erg s . The dynamics is identical 
to the iMcCrav fc Kafatosl (|1987|) model and hence can be inverted to give, 

t 7 = (i? s /267pc)(^ s /15.7km s^ 1 )" 1 (A3) 

and 

(Wno) = (^s/267 pc) 2 (F s /15.7 km s" 1 ) 3 , (A4) 

which provide the required input parameters as a function of the observed properties of the supershell. The expression 
for the age is identical to that derived earlier in this work. Substituting the values for Vs and no, the required 
mechanical luminosity for the cluster is 3.6 x 10 38 erg s _1 . 

This mechanical luminosity may now be compare d directly with those obtained from population synthesis models. 
For the fiducial Starburst99 (|Leitherer et al.|[l999l) simulation described earlier in this work, the mean mechanical 
luminosity in the first 15 Myr from stellar winds and supernovae is given as 2.57 x 10 34 erg s _1 Mq . Comparing this 
with the energy budget, requires a massive ~ 1.4 X 10 M Q young stellar cluster. But, this is lower than the mass 
required in the earlier analysis and hence the observed UV flux is consistent with a cluster massive enough to drive 
the supershell. 



